UMN-TH-3105/12, FTPI-MINN-12/19 

The variation of fundamental constants and the role of A = 5 and A = 8 nuclei on 

primordial nucleosynthesis 

Alain CocB 

Centre de Spectrometrie Nucleaire et de Spectrometrie de Masse (CSNSM), 
IN2P3-CNRS and Universite Paris Sud 11, UMR 8609, Bat. 104, 91405 Orsay Campus (France) 

Pierre DescouvemontQ 

Physique Nucleaire Theorique et Physique Mathematique, CP. 229, 
Universite Libre de Bruxelles (ULB), B-1050 Brussels, Belgium 

Keith A. OliveS 

William I. Fine Theoretical Physics Institute, School of Physics and Astronomy, 
University of Minnesota, Minneapolis, MN 55455 (USA) 

Jean-Philippe Uzard and Elisabeth Vangion^] 

Institut d'Astrophysique de Paris, UMR-1095 du CNRS, 
Universite Pierre et Marie Curie, 98 bis bd Arago, 75014 Paris (France) 

We investigate the effect of a variation of fundamental constants on primordial element production 
in big bang nucleosynthesis (BBN). We focus on the effect of a possible change in the nucleon-nucleon 
interaction on nuclear reaction rates involving the A — 5 ( 5 Li and 5 He) and A — 8 ( 8 Be) unstable 
nuclei and complement earlier work on its effect on the binding energy of deuterium. The reaction 
rates for 3 He(d,p) 4 He and 3 H(d,n) 4 He are dominated by the properties of broad analog resonances 
in 5 He and 5 Li compound nuclei respectively. While the triple alpha process 4 He(aa, j) 12 C is 
normally not effective in BBN, its rate is very sensitive to the position of the "Hoyle state" and 
could in principle be drastically affected if 8 Be were stable during BBN. The nuclear properties 
(resonance energies in 5 He and J Li nuclei, and the binding energies of 8 Be and D) are all computed 
in a consistent way using a microscopic cluster model. The n(p,7)d, 3 He(d,p) 4 He, 3 H(d,n) 4 He 
and 4 He(aa, j) 12 C, reaction rates are subsequently calculated as a function of the nucleon-nucleon 
interaction that can be related to the fundamental constants. We found that the effect of the 
variation of constants on the 3 He(d,p) 4 He, 3 H(d,n) 4 He and 4 He(aa,7) 12 C reaction rates is not 
sufficient to induce a significant effect on BBN, even if 8 Be was stable. In particular, no significant 
production of carbon by the triple alpha reaction is found when compared to standard BBN. We 
also update our previous analysis on the effect of a variation of constants on the n(p,7)d reaction 
rate. 
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I. INTRODUCTION 

Constraints on the possible variation of fundamental 
constants are an efficient method of testing the equiva- 
lence principle [l|, [2] , which underpins metric theories of 
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gravity and in particular general relativity. These con- 
straints are derived from a wide variety of physical sys- 
tems and span a large spectrum of rcdshifts and physi- 
cal conditions, from the comparison of atomic clocks in 
the laboratory, the Oklo phenomena, to quasar absorp- 
tion spectra up to a typical redshift of order z ~ 2 and 
big bang nucleosynthesis (BBN) at a redshift of order 
z ~ 10 9 . 

Primordial nucleosynthesis is considered a major pillar 
of the standard cosmological model (see e.g., Refs. @n. 
Using inputs from WMAP for the baryon density j4|, 
BBN yields excellent agreement between the theoreti- 
cal predictions and astrophysical determinations for the 
abundances of D and 4 He despite the discrepancy 

between the theoretical prediction of 7 Li and its deter- 
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mined abundance in halo stars [8|. Indeed, BBN has 
been used extensively to constrain deviations from the 
standard model framework, and in particular from gen- 
eral relativity, see e.g., Ref. [Io| . 

The effects of the variation of fundamental constants 
on BBN predictions is difficult to model because of the 
intricate structure of QCD and its role in low energy nu- 
clear reactions and because one cannot restrict the anal- 
ysis to a single constant. One can, however, proceed in a 
two step approach: first by determining the dependencies 
of the light element abundances on the BBN parameters 
and then by relating those parameters to the fundamen- 
tal constants (see Section 3.8 of Ref. [l[ for an up-to-date 
overview). While early works have mostly focused on a 
single parameter such as the fine structure constant fill ] , 
the Higgs vacuum expectation value (vev) [l2l - fl4j or the 
QCD scale [lj| , in many theories which allow for the vari- 
ation of fundamental parameters, often, the variation of 
sever al p arameters are correlated in a model dependent 
way 

The variation of a fundamental parameter such as 
the fine structure constant will affect the BBN analy- 
sis through the proton-to-neutron mass difference and 
the neutron lifetime (ll| as well as the deuterium bind- 
ing energy [l8| and the binding energies of other light 
nuclei such as tritium, helium-3 and 4, lithium-6 and 7, 
and beryllium-7 [l9j]. These effects can in principle be 
used to probe the coupled variations of several parame- 
ters [IchII. 

Following our previous work (22| - [24| . we allow for a 
variation of all fundamental constants and in order to re- 
duce the arbitrariness, we focus on scenarios in which the 
variations of the different constants are correlated. In ef- 
fectively all unification models of non-gravitational inter- 
actions, a variation in the fine structure constant is asso- 
ciated with the variation in other gauge couplings PUGS- 
Any variation of the strong gauge coupling a s will induce 
a variation in the QCD scale, Aqcd, as can be seen from 
the expression 
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valid for a renormalization scale fi > nit- In this ex- 
pression m Ci b,t are the masses of the charm, bottom and 
top quarks. Since the masses of the quarks are propor- 
tional to the product, hv, of a Yukawa coupling h and 
Higgs vacuum expectation value u, any variation of the 
Yukawa couplings will also induce a variation of Aqcd- 
These variations can be related by 
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The coefficient R is determined by the particular grand 
unified theory and particle content of the theory which 
control both the value of «(Mgut) = a s (A/gut) and 
the low energy relation between a and a s , leadin g t o 
a considerable model dependence in its value (2|| |26| . 



Here we shall assume a typical value, R ~ 36 [l7j, l27| . 
Furthermore, in theories in which the electroweak scale 
is derived by dimensional transmutation, changes in the 
Yukawa couplings (particularly the top Yukawa, h t ) lead 
to exponentially large changes in the Higgs vev. In such 
cases, the Higgs expectation value is related to the Planck 
mass, Mp, by 
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where c is a constant of order unity, and at = h^/4ir. 
Thus we can write, 
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and, as in Ref. [22[, we take S ~ 240, though there is 
considerable model-dependence in this value as well. For 
example, in supersymmetric models, S can be related 
to the sensitivity of the Z gauge boson mass to the top 
Yukawa, and may take values anywhere from about 80 to 
500 [28| . This dependence gets translated into a variation 
in all low energy particle masses [l2| . In addition, in 
many string theories, all gauge and Yukawa couplings 
are determined by the expectation value of a dilaton. 
Therefore, once we allow a to vary, virtually all masses 
and couplings are expected to vary as well, typically much 
more strongly than the variation induced by the Coulomb 
interaction alone. 

The use of coupled variations has led to significantly 
improved constraints in a wide range of environments 
ranging from big ban g n ucleosynthesis p7l -l 2C j . I22I . I29I - 
[32( 1 , the Oklo reactor l33l. meteoritic data [2l|, l33l . the 
microwave background |3lL IHH , stellar evolution |23| and 
atomic clocks [24j. Concerning BBN, the effect of cou- 
pled variations has mostly focused on the binding energy, 
Bjj, of deuterium [T"' 
lated investigations) 
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[22|, |29| (see also Ref. 
The importance of Bu is easily 
understood by the fact that the equilibrium abundance 
of deuterium and the reaction rate p{n, 7)D both depend 
exponentially on Bd and on the fact that deuterium is in 
a shallow bound state. Indeed, in Ref. (22|, we found that 
even a relatively small variation in the gauge or Yukawa 
couplings of order of a few xlO -5 had a significant effect 
on the light element abundances. In particular, using [22[ 
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a variation in the Yukawa couplings of 2 x 10~ 5 induces 
a relative variation in Br, of about 4%. By decreasing 
Bo, nucleosynthesis begins later at a lower temperature 
ultimately suppressing the 7 Li abundance. 

It is well known in principle, that the mass gaps at 
A =5 and A = 8, prevent the nucleosynthetic chain from 
extending beyond 4 He. Although some 6 Li and 7 Li is 
produced, their abundances remain far below that of the 
lighter elements, while B, Be, and CNO isotopes are pro- 
duced in even smaller amounts. The presence of these 
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gaps is caused by the instability of 5 He, 5 Li and 8 Be with 
respect to particle emission: their lifetimes are as low as a 
few 1CT 22 s for 5 Hc and 5 Li and « 10~ 16 s for 8 Bc. More 
precisely, 5 He, 5 Li and 8 Be are respectively unbound by 
0.798, 1.69 and 0.092 MeV with respect to neutron, pro- 
ton and a particle emission 1 . Variations of constants will 
affect the energy levels of the unbound 5 He, 5 Li and 8 Bc 
nuclei [IH, H(| and hence, the resonance energies whose 
contributions dominate the reaction rates. In addition, 
since 8 Be is only slightly unbound, one can expect that 
for even a small change in the nuclear potential, it could 
become bound and may thus severely impact the results 
of standard BBN (SBBN), in a similar way that a bound 
dineutron impacts BBN abundances [IB] • It has been sus- 
pected that stable 8 Be would trigger the production of 
heavy elements in BBN, in particular that there would 
be significant leakage of the nucleosynthetic chain into 
carbon. Indeed, as we have seen previously j23j . changes 
in the nuclear potential strongly affects the triple alpha 
process and as a result, strongly affects the nuclear abun- 
dances in stars. 

This article investigates in detail the effect of the vari- 
ation of fundamental constants on the properties of the 
compound nuclei 5 He, 5 Li, 8 Be and 12 C involved in the 
3 H(d,n) 4 He, 3 He(d,p) 4 He and 4 Hc(aa, 7 ) 12 C reactions, 
and their consequences on reaction rates and BBN abun- 
dances. In addition, we consider the particular case of 
stable 8 Be. 

In Section UH we briefly present the microscopic clus- 
ter model used to determine resonance properties. Sec- 
tion IIIII focuses on 8 Be and on the CNO production 
for the cases of both unbound and bound 8 Bc. We 
compute the C/H ratio as a function of the parameter 
5nn- Section HVl focuses on the 4 He production by the 
3 Hc(d,p) 4 Hc and 3 H(d,n) 4 Hc reactions. Each of these re- 
actions contains a broad s-wave resonance at low energies, 
and their reaction rates may depend on the resonance en- 
ergies. Section IVl summarizes the BBN constraints on the 
variation of the nuclear interaction, hence extending our 
previous analysis [22j ■ Section I VII provides a summary 
and our conclusions. 



II. OUTLINE OF THE NUCLEAR MODEL 

We follow the formalism introduced in Ref. [23| in or- 
der to model the effect of the variation of the nucleon- 
nucleon (N-N) interaction. We adopt a phenomenologi- 
cal description of the different nuclei based on a cluster 
model in which the wave functions are approximated by 



1 Although 8 Li and 8 B also contribute to the mass gap due to 
their short lifetimes (on BBN timescales), they are more deeply 
unbound, and a far greater change in the fundamental couplings 
would be needed to affected their stability. We will not consider 
them further here. 



clusters of two or three a wave functions. In a micro- 
scopic theory, the wave function of a nucleus with nu- 
cleon number A, spin J, and total parity 7r is a solution 
of a Schrodingcr equation with a Hamiltonian given by 

A A 

i— 1 z>j — 1 

where Tj is the kinetic energy of nucleon i, and Vij a 
nucleon-nucleon (N-N) interaction. In general, the po- 
tential depends on space, spin and isospin coordinates of 
nucleons i and j, and can be decomposed as 

Vij = Vg + Vg, (2.2) 

where V c and V N represent the Coulomb and nuclear 
interactions, respectively. 

Solving the Schrodinger equation associated with 
Hamiltonian (|2.1[) is a difficult problem, in particular for 
nuclear reactions. We use a microscopic cluster model, 
where the nucleons are assumed to form " groups" , called 
clusters. This approximation is known as the Resonating 
Group Method (RGM) [37j], and the wave function of a 
two-cluster system is approximated as 

qjMn = A^ 2 g J ' r {p)Yf(fl), (2.3) 

where $i and $2 are the internal wave functions of the 
clusters (defined in the shell model), and A is the A- 
nucleon antisymmctrizor, which accounts for the Pauli 
principle. In Eq. (|2.3p . g' ,7T (p) is the relative function, 
depending on the relative coordinate p, and to be deter- 
mined from the Schrodingcr equation. With these wave 
functions, any physical quantity, such as spectroscopic 
properties of the nucleus, or nucleus-nucleus cross sec- 
tions, can be computed. 

Eq. (|2.3I) is written for two clusters, but the ex- 
tension to three clusters is feasible (38|. Currently, 
more sophisticated microscopic models (such as the No 
Core Shell Model I39l| or the variational Monte Carlo 
(VMC) method [40() are available for few- nucleon sys- 
tems. These models use realistic interactions (such as 
Argonne AV18 [4l[) but are quite difficult to apply with 
nucleus-nucleus collisions. In contrast, the RGM is well 
adapted to the nuclear spectroscopy and to reactions, 
but the use of simple cluster wave functions (such as 
the a particle which is described by four 0s orbitals) 
makes it necessary to adopt effective N-N interactions. 
We use here the Minnesota potential (42[, wen adapted 
to low-mass systems. This central potential reproduces 
the experimental dcuteron energy. It simulates the miss- 
ing tensor force by an appropriate choice of the central 
interaction. The Minnesota interaction V^v is described 
in detail in Ref. H|. 

To take into account the variation of the fundamen- 
tal constants, we introduce the parameters 5 a and 5nn 
to characterize the change of the strength of the electro- 
magnetic and nucleon-nucleon interactions respectively. 
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This is implemented by modifying the interaction poten- 
tial flUD so that 



III. PRIMORDIAL CNO PRODUCTION AND 

8 BE 



(i + & a )vS + (i + s NN )Vi 
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(2.4) 



Such a modification will affect Bp, the energy levels 
of 8 Be and 12 C simultaneously, as well as the resonant 
reactions involving A = 5 nuclei such as 3 He(d,p) 4 He 
and 3 H(d,n) 4 Hc. If 8 Be becomes bound, one will need 
to calculate the two reaction rates 4 He(a,7) 8 Be and 
8 Be(a, 7 ) 12 C. 

In Ref. [H, we investigated the 8 Be and 12 C(0^) en- 
ergies by scaling the Minnesota interaction. The Snn 
parameter then provides a link, for the Minnesota poten- 
tial, between the deutcron energy Bd and the 2a and 
3a resonance energies. Here we extend this idea to the 
3 H(d,n) 4 He and 3 He(d,p) 4 He reactions. However both 

reactions are known to be dominated by a low energy ^ 
resonance (at E e ^ p = 0.048 MeV for 5 He and E e * p = 0.21 
MeV for 5 Li). From simple angular-momentum couplings 
it is easy to see that this resonance corresponds to an s 
wave in the entrance channel, and to a d wave in the exit 
channel. Consequently the coupling between these chan- 
nels can be described by a tensor force only. As men- 
tioned earlier, this component is neglected in the Min- 
nesota interaction. However, as our main interest is the 
variation of the resonance energy, we used single channel 
3 He+d and 3 H+d approximations, and performed vari- 
ous calculations by modifying Snn (see Section III). This 
approximation is justified by the fact that these reso- 
nances essentially have a 3+2 structure. 

Snn is a phenomenological parameter that can be re- 
lated to the fundamental constants through the depen- 
dence of the deuterium binding energy on Snn- With 
the Minnesota N-N interaction, we find 



AB D /B D = 5.716 x 5 



NN ■ 



(2.5) 



Note that other forces may provide a slightly different 
dependence. However, a consistent treatment of Bu 
and of resonance properties in 5 Li and 5 He requires the 
same effective interaction, such as the Minnesota po- 
tential. Thus we have the possibility of relating Snn 
to the gauge and Yukawa couplings if one matches this 
prediction to a potential model via the a and u> meson 
masses fTpL I22I I2II [43| or the pion mass, as suggested in 
Refs. Il4 144 1451 . In Ref. [Hf, it was concluded that 
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which led to the expression in Eq. (| 1 . 5[) . Eq. (|2.5|) 
can then link any constraint on Snn to the three 
fundamental constants (h s ,v,N). 



CNO production in SBBN has been investigated in 
Ref. and most recently revisited in Coc et al. [§]. 
The direct detection of primordial CNO isotopes seems 
highly unlikely with the present observational techniques 
but it is important for other applications. In particu- 
lar, it may significantly affect the dynamics of population 
III (Pop. Ill) stars since hydrogen burning in low mass 
Pop. Ill stars proceeds through the slow pp chains until 
enough carbon is produced, through the triple alpha re- 
action, to activate the CNO cycle. The minimum value of 
the initial CNO mass fraction that would affect Pop. Ill 
stellar evolution was estimated to be 10 -10 [46| or even 
as low as 10~ 12 for less massive stars [47|- This is only 
two orders of magnitude above the SBBN CNO yields 
obtained using current nuclear reaction rates. The main 
difficulty in BBN calculations up to CNO is the extensive 
network (more than 400 reactions) needed, including n, 
p, a, but also d, t and 3 Hc, induced reactions on both 
stable and radioactive targets. 

CNO production (mostly 12 C) in SBBN is found to be 
in the range CNO/H = (0.2 - 3.) x 10~ 15 % In sce- 
narios with varying constants, this number needs to be 
compared with the 4 Hc(aa, 7) 12 C production of 12 C, in 
particular if its rate is dramatically increased by a vari- 
ation of the 8 Be ground state and Hoyle state position 
as considered by Ekstrom et al. [23|. Since it is a res- 
onant reaction, its rate is very sensitive to the nuclear 
interaction and we recall that a 1.5% variation of the N- 
N interaction (—0.009 < Snn < +0.006) would induce a 
change in the rate between «18 to sw 2 orders of magni- 
tude for temperatures between T = 0.1 to 1.0 GK (see 
Fig. 3 in Ref. (H). 

In this section, we investigate the effect of 8 Be on 
primordial CNO production. We consider two cases in 
which 8 Be is either unbound (§ IIII Ap or bound (§ 1111 B[) . 
We then derive the BBN predictions (§ IIII C[) in each 



case. 



A. Unbound 8 Be 

When the N-N interaction is modified by less than 
0.75% (i.e. Snn < 7.52 x 10~ 3 ), 8 Be remains unbound 
w.r.t. two a-particle emission. We can therefore take the 
4 He(aa, 7) 12 C rate as a function of Snn as calculated by 
Ekstrom et al. [23| . We recall that, in the framework 
of the cluster model using the Minnesota interaction, we 
obtained that the energy of the 8 Be ground state with 
respect to the a + a threshold was given by [23j 

-B 8 = E R ( s Be) 

= (0.09184- 12.208 x S NN ) MeV (3.1) 

where Bg, is the 8 Be binding energy with respect to two 
alpha break-up (with this convention, £>g < for un- 
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bound 8 Bc). We recall that, within the same model, we 
obtained 

E R { 12 C) = (0.2876 - 20.412 x 5 NN ) MeV (3.2) 
for the dependence of the Hoyle state resonance. 

B. Bound 8 Be 

When Snn ^ 7.52 x 10~ 3 , 8 Be becomes bound 
and should be considered as a stable isotope during 
BBN. Hence, we have to calculate two reaction rates: 
4 He(a,7) 8 Be and 8 Be(a,7) 12 C. The calculation of the 
rate of the second reaction can be achieved using the 
sharp resonance formula [48j with the var ying parame- 
ters of the Hoyle state from Ekstrdm et al. [23|. For the 
first reaction, 4 He(ct,7) 8 Be, we have performed a dedi- 
cated calculation using the potential by Buck et al. [49| 
to obtain the astrophysical S- factor displayed in Figure Q] 
for values of the 8 Be binding energy of B s = 10, 50 and 
100 kcV. The Buck potential is expressed as a Gaussian, 
and accurately reproduces the experimental a + a phase 
shifts up to 20 MeV. The initial 2+ and final 0+ wave 
functions arc computed in the potential model, assuming 
that the ground state is slightly bound. The cross sec- 
tion is then determined from integrals involving the wave 
functions and the E2 operator (see Ref. [48| for detail). 
The broad structure corresponds to the well known 2 + 
resonance in a — a scattering and it is easily concluded 
from Figure [1] that the S- factor remains relatively insen- 
sitive to a change in 5nn- It was shown in Ref. [50J that 
describing 8 Be as a bound state (Er < 0) or as a low- 
energy resonance (Er > 0), has a small effect on the 
a(a, 7) 8 Be capture cross section. 

Figure [2] depicts the reaction rate for B% = 10 and 100 
keV relative to the case with B$ = 50. The rate depends 
very little on the 8 Be binding energy for £?§ > and 
the rate changes by less than ~10% for the three values 
of B$ considered. As a result, we can safely neglect the 
difference in the rates once B$ > 0. The reaction rate is 
essentially given by the radiative capture cross section at 
the Gamow energy Eq(T), which is proportional to E^, 
where the photon energy is E 1 = E cm + B 8 and E cm is 
the alpha-alpha center-of-mass energy. 
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FIG. 1: The astrophysical S-factor for the 4 He(«,7) 8 Be 
reaction, assuming that 8 Be is bound by 10 (blue), 50 
(red) and 100 (green) keV, corresponding to 5 NN = 
0.0083, 0.0116, 0.0156 respectively, from bottom to top. 
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FIG. 2: The relative variation of the 4 He(a,7) 8 Be reaction 
rate assuming that 8 Be is bound by 10 (blue), 50 (red) and 
(green) 100 keV, relative to the 50 keV rate. 



C. BBN calculations 

CNO is produced at a very low level in SBBN. The 

chain leading to carbon is dominated by the following 
reactions: 

7 Li(a, 7 ) n B 7 Li(n,7) 8 Li(a,n) n B (3.3) 
followed by 

n B(p, 7 ) 12 C n B(d,n) 12 C, 

n B(d,p) 12 B n B(n, 7 ) 12 B (3.4) 



which bridge the gap between the A < 7 and A > 12 
nuclei @, M- 

To disentangle standard CNO production through the 
reactions listed above with the one proceeding through 
the triple-alpha reaction, we reduced the network to 
the 15 reactions involved in A < 8 nucleosynthesis plus 
4 He(aa, 7) 12 C, i.e. we turned off the other reactions, in 
particular those listed in Eqs. (|3.3[) and (|3.4[) . As we 
are investigating a possible enhancement of CNO pro- 
duction we only considered positive Snn values that lead 
to a higher triple-alpha reaction rate. 

The CNO yield as a function of 5nn is displayed in 
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Figure [3l The carbon abundance shows a maximum 
at Snn ~ 0.006, C/Hw 10 -21 , which is six orders of 
magnitude below the carbon abundance in SBBN. This 
can be understood as follows: Figure 2] displays the 
4 Hc(aa, 7) 12 C rate as a function of Snn for tempera- 
tures relevant to BBN, i.e., from 0.1 to 1 GK. As one 
can see, the variation of the rate with Snn is limited 
at the highest temperatures where BBN production oc- 
curs so that the amplification of 12 C production does 
not exceed a few orders of magnitudes. Indeed, while 
stars can process CNO at 0.1 GK over billions of years, in 
BBN the optimal temperature range for producing CNO 
is passed through in a matter of minutes. This is not 
sufficient for 12 C (CNO) nucleosynthesis in BBN. Fur- 
thermore, the baryon density during BBN remains in the 
range 10~ 5 to 0.1 g/cm 3 between 1.0 and 0.1 GK, sub- 
stantially lower than in stars (e.g. 30 to 3000 g/cm 3 in 
Pop. Ill stars). This makes three-body reactions like 
4 He(aa, 7) 12 C much less efficient compared to two-body 
reactions. Finally, in stars, 4 He(aa, 7) 12 C operates dur- 
ing the helium burning phase without significant sources 
of 7 Li, d, p and n that allow the processes listed in Eqs. 
(USD and dm 
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0.001 0.002 0.003 0.004 0.005 0.006 0.007 

S NN 

FIG. 4: The 4 He(aa, 7) 12 C rate as a function of 5jvjv at 
constant temperature relevant for BBN from 0.1 GK (lower 
curve) to 1 GK (upper curve), by 0.1 GK steps. 
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FIG. 3: 12 C production (in number of atoms relative to H) 
through the 4 He(aa, 7) 12 C reaction (left) or via a stable 8 Be 
(right) as a function of the N-N interaction. For clarity, the 
rates of all other CNO producing reactions are set to zero to 
study these specific channels. The dotted line just connects 
the results of the two types of calculations: via 4 He(aa, 7) 12 C 
a s in |p and S fllL&l or 

§inm 



1 He(a,7) 8 Be staWe (a, 7 ) 12 C as in 



channel is given by 



\ av ' <* "F7W\ exp 
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k B T 



(3.5) 



r 7 the exit 



where r Q is the entrance (alpha-)width 
(gamma-)width and T is the total width (r = T a + if 
these are the only open channels). Figure 1 in Ref. (23| 
displays the 12 C level scheme: the radiative width, r 7 , is 
associated with the decay to the first 12 C excited state 
at 4.44 MeV as the decay to the ground state proceeds 
only through the much less efficient electron-positron 
pair emission. The corresponding decay energy is then 
Ey(E R ) = 3.21 MeV + AE R { 12 C). Equation E3] shows 
that for a fixed Er, i.e. Snn, the contribution increases 
with temperature as seen in Figs. 2] and [5] While the 



radiative width T 1 (Er) 



is almost insensitive 



to Er(6nn), ^oi(Er) is very sensitive to Er(Snn) vari- 
ations because of the Coulomb and centrifugal barriers 
penetrability, Pg(E). The reduced widths 7 2 , defined by: 



T X (E) = 2 1 2 X P ( (E) (^ 7 ), 



(3.6) 



The maximum of the 12 C production as a function of 
6nn in Figure [3] reflects the maxima in the 8 Be(a, 7 ) 12 C 
and 4 He(aa, 7 ) 12 C rates displayed in Figs. 2] and [SJ They 
are due to the effect of the sharp resonances (in both 
the 8 Be ground state and the 12 C "Hoyle state" that 
dominate the cross section). 

The contribution of a sharp resonance in the (a, 7 ) 



are corrected for these effects so that they reflect the nu- 
clear properties only, and are, as a good approximation, 
independent of Snn- 

Depending on whether r Q -C T 7 or T 7 -c T a , the 
sensitivity of (av) (Eq. 13. 5[) to Er or Snn variations is 
very different. This is due to the very different en ergy 
dependence of T a and T 7 , as discussed in detail in |5l| . 
In the latter case, using Eq. (|3.5[) . the sensitivity of the 



7 



rate to Er (5 



n n ) variations is simply given by : 

dln(av) Er 
d\nE R = ~k~~f 

is reduced to I\ 



(3.7) 

as the prefactor in Eq. (|3.5[) is reduced to T 7 which is 
almost constant. Since SEr and Snn have opposite signs 
(Eqs. 13.11 and 13. 2|) . the rate increases with Snn- In the 
former case, the same factor is reduced to the very energy 
dependent r Q and we have: 



(<rv) oc 7 2 exp 




Er k B T 



(3.8) 



where the penetrability, Pi(E), has been approximated 
by exp(— y^Ec/E) with Gamow energy, Eq. It is well 
known that the exponential in Eq. (|3.8[) can be well ap- 
proximated (see e.g. Ref. [48[) by 



exp 



Er — Eq 
AEo/2 



(3.9) 



with 



/i\V3 fne 2 Z 1 Z 2 kT 



2/3 



27 V ft 
0.1220 (Z?Zf A) 1/3 T 9 2/3 MeV (3.10) 



and 



AE = 4 (£ kT/3) 1/2 

= 0.2368 (Z 2 Z 2 A) 1/6 T 9 5/6 MeV. (3.11) 

that define the Gamow window. Recalling that the re- 
duced width 7 2 only reflects the nuclear structure and is 
assumed to be constant, it is straightforward to calculate 
the sensitivity of the rate to Er (Snn) variations: 



d\n(av) 
d\nE R 



E (T) - E R (6 NN ) 
AE (T)/2 



(3.12) 



Since, for large Snn and T, we have Eq > Er, the rate 
decreases with Snn- 

The condition T a = T 7 that marks the boundary be- 



tween these two opposite evolutions in Eq. fl3.5[) can be 
found in Fig. A.l of Ekstrom et al. [2j|, at Snn ~ 0.006. 
It also corresponds to the maximum of the 8 Be(a,7) 12 C 
rate depicted in Fig. [5J The more complicated two step, 
three body 4 He(aa,7) 12 C reaction shows a similar de- 
pendence. 

To summarize, for 5nn ~ 0.006, the rate decreases 
(increases) as a function of Er (Snn) because of the 
dominating exponential factor, exp(— Er/YT), while for 
Snn > 0.006, it increases (decreases) because of the pen- 
etrability. This evolution is followed by the 12 C produc- 
tion displayed in Fig. [3] 

For Snn > 0.00752, when 8 Be is bound, 12 C produc- 
tion drops to C/H w 5 x 10~ 23 for B$ = 10. For larger B$, 



the abundance drops sharply as seen in Figure [3] For B% 
= 50 keV, C/H ~ 5 x 10~ 29 and is no longer in the range 
shown in the figure. For B$ = 100 keV, corresponding to 
Snn = 0.0156, the Hoyle state is even below threshold 
and the production is vanishingly small. If 8 Be is bound, 
reactions that normally produce two a-particles could 
form 8 Be instead. We considered the following reactions 



7 Bc(n,7)2a 
7 Li(n, 7 ) 8 Li(/? + )2a 
7 U(d,n)2a 
7 Be( 3 He, 2p)2a 
7 Li( 3 He, d)2a 
7 Li(i,2n)2a 



7 Li(p, 7 )2a 
7 Be(d,p)2a 
7 Bc(i, np)2a 
7 Be(n,7)2a 
7 Be(t,d)2a 
7 Li( 3 He,np)2a 



using the same rates as in Ref. @ but replacing 2a by 
8 Be. The only significant enhancement comes from the 
7 Li(d,n)2a reaction but even in the most favorable case 
(Bg = 10 keV), C/H reaches w 10" 21 . This is still six 
orders of magnitude below the SBBN yield Q that pro- 
ceeds via the reactions listed in Eqs. f|3.3[) and (|3.4p . 



T = 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, 1.0 GK 




0.002 0.004 0.006 0.008 



0.012 0.014 



FIG. 5: As in Fig. g] but for the s Be(a,7) 12 C rate. Above 
Snn ~ 0.006, the contribution of the Hoyle state decreases 
(see text), leaving only the contribution of the 3~ resonance 
at Er = 2.274 MeV at the highest Snn values. 

Figure |JJ] displays the evolution of the 12 C and 8 Be 
mass fractions as a function of time. They both increase 
with time until the 4 Hc(a,7) 8 Be drops with decreas- 
ing temperature. Afterwards, equilibrium between two 
a-particle fusion and 8 Be photodissociation prevails as 
shown by the dotted lines. Indeed, the Be mass fraction 
is 



2 R 



oc 



Ti exp (-B 8 /k B T) 



(3.13) 
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where we have taken the reverse ratio, R, to be pro- 
portional to exp (Bg/ksT 1 ) xT? and the baryon density, 
p oc T 3 . For the highest values of B$, the 8 Be mass frac- 
tion increases until, due to the expansion, equilibrium 
drops out, as shown by the late time behaviour of the 
upper curve (B$ = 100 keV) in Figure [6] 



Be 10, 50 and 100 keV bound 

io 3 io 4 




Time (s) 

FIG. 6: 12 C and 8 Be mass fractions as a function of time, 
assuming 8 Be is bound by 100, 50 and 10 keV as shown by 
the upper to lower curves respectively. (Only the 12 C mass 
fraction curve, for _Bg = 10 keV, is shown; others are far 
below the scale shown). The dotted lines correspond to the 
computation at thermal equilibrium. 

In conclusion, if one keeps only the first 15 nuclear 
reactions while adding the possibility of primordial car- 
bon production through 3a-process, one gets a typical 
abundance of only order C/H^ 10 -21 which remains six 
orders of magnitude smaller than the SBBN carbon pro- 
duction, which is of order C/H^ 10 -15 . This explicitly 
shows that CNO production in BBN cannot be increased 
through cither an increase of the 3ev-reaction rate or even 
through the stabilization of 8 Be. 

IV. REACTIONS INVOLVING A = 5 NUCLEI 

A. Introduction 

In contrast to the CNO elements, the 4 He, D, 3 He 
and 7 Li isotopes are produced in observable quantities 
in BBN and can thus constrain <5jvjv- It would be desir- 
able to have the dependence of each of the main SBBN 



reactions to the N-N interaction or other fundamental 

quantities. This was achieved in Ref. [22j for the first 

two BBN reactions: the n-o-p weak interaction and the 

p(n,7)d bottleneck. 

Here, we propose to extend this analysis to the 

3 H(d,n) 4 Hc and 3 He(d,p) 4 He reactions that proceed 

through the A = 5 compound nuclei 5 He and 5 Li. In 

these two reactions, the rates are dominated by the can- 
't + 

tribution of a | resonance whose properties can be cal- 
culated within the same microscopic model that we used 
for 4 Hc(aQ!, 7) 12 C, but with 3 H+d and 3 He+d cluster 
structures. Unlike the case for 8 Be, the lifetime of the 
5 He and 5 Li states is extremely short (the width of the 
8 Be ground state is 6 eV, whereas the widths of the § 
resonances in 5 He and 5 Li are of the order of 1 McV). 
Therefore the issue of producing A = 5 bound states, 
or even shifting their ground state energy down to the 
Gamow window, is not relevant. Even a two step pro- 
cess, like the triple alpha reaction, where 5 He or 5 Li in 
thermal equilibrium would capture a subsequent nucleon 
to form 6 Li is completely negligible because they are un- 
bound by ~1 MeV compared with the 92 keV of 8 Bc. 
Hence no significant equilibrium abundance of A = 5 nu- 
clei can be reached. 

The analysis of the effect of 5mn on the other rates 
requires additional effort, but should be smaller because 
important reactions like 3 Hc(a, 7) 7 Bc do not display res- 
onant behaviour and the S'-factor does not change much 
with energy (see e.g. Fig. Hi in Ref. [12 )■ In order to 
proceed, we used the following procedure 

• We performed single-channel 3 H+d and 3 He+d cal- 
culations of the resonance energy as a function of 
5nn- 

• We analyzed the experimental cross sections within 
an i?-matrix approach. 

• With the energy dependence obtained from the first 
step, we computed the reaction rate for various Snn 
or Bu values. 



B. Resonance energy in 5 He and 5 Li 

In principle an RGM calculation of the 3 H(d,n) 4 Hc and 
3 Hc(d,p) 4 He cross sections could be performed as it was 
done for the 2a and 3a systems. The main difference is 
that each reaction involves two channels, and their cross 
section is dominated by a low-energy | + resonance. As 
mentioned before, a central N-N interaction, such as the 
Minnesota potential, does not provide a coupling between 
the 3 H+d and a+n (or mirror) channels, since the ten- 
sor force is missing. Consequently, the neutron width 
is strictly zero, and obtaining a realistic S'-factor with 
the Minnesota interaction is not possible. We therefore 
use the RGM to estimate the sensitivity of the resonance 
energy as a function of 6 n n- 
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Using the parameterization (|2.4|) for the nucleon- 
nucleon interaction, we modify the resonance energy. 
Both the excitation energies of the | + resonance and of 
the thresholds vary. We find 



AE R = -0.327 x 8 



NN 



for 3 H(d,n) 4 He and 



-0.453 x 8 NN 



(4.1) 



(4.2) 



for 3 He(d,p) 4 He (units are MeV). These energy depen- 
dences are much weaker (~ 20-30 keV for \Snn\ < 0.03) 
than for 8 Bc and 12 C (sec Eqs. fl3H]) and dHlj)). This is 
expected for broad resonances which are weakly sensitive 
to the nuclear interaction [37| . 

In contrast, Berengut et al. (30j find a stronger energy 
dependence. These authors perform VMC calculations 
with realistic N-N interactions, which provide better d 
and 3 H/ 3 He wave functions. However the VMC approach 
is not well adapted to broad resonances, such as those 
observed in 5 He and 5 Li. Our cluster model, although 
using a simpler N-N interaction, is more suited to un- 
bound states since the asymptotic Coulomb behaviour of 
scattering states is exactly taken into account. 



i?-matrix fits of the 3 H(d,n) 4 He and 
3 He(d,p) 4 He S-factors 



To be consistent with our previous work, we want to 
reproduce, for 8nn=Q, the ex per imental 5-factors ob- 
tained (see references in Ref. j52|) by a full i?-matrix 
analysis, but for convenience we restrict ourselves here 
to the single pole i?-matrix approximation which will be 
shown to be sufficient. For the 3 H(d,n) 4 He reaction, we 
use the parameterization of Barker [53J , which reproduces 
the resonance corresponding to the corresponding to the 

r, _l_ 

| state at 16.84 MeV. Because the widths are energy 
dependent, the maximum of the cross section at E R P = 
0.048 MeV, differs from the minimum of the denominator 
that appears in the single pole i?-matrix prescription: 



a{E) oc 



{hef 



r in (£)r out (£) 



liE {E* R + AE R -E) 2 +T 2 (E)/4 



(4.3) 



where /i the reduced mass, T = T ln + T out and we 
have dropped numerical factors. The astroph ysical S- 
factor displayed in Figure [7] is just a(E)Eexp \J (Eq/E) 
where the exponential factor approximates the en- 
ergy dependence due to the penetrability. The widths, 
r x , are related to the reduced widths 7^ by Eq. (|3.6p . 
and AE%,(Er) is given by Eqs. (5.5) in Descouvcmont 
& Baye (5J]- Note that AE R is the shift factor of R- 
matrix theory, not related to the variation of constants, 
and hence is different from the AEr(Snn), that we will 
discuss below. Using E R = 0.091 MeV, the reduced en- 
trance width 7^ = 2.93 MeV and the reduced exit width 
7^ = 0.0794 MeV provided by Barker [53|, our results 



(no t a fit) are in perfect agreement with the i?-matrix fit 
[52j of the experimental data shown in Figure [7] that we 
used in previous work 0, [t| [22|, [55| . 

Barker was also able to reproduce the existing experi- 
mental data for the 3 Hc(d,p) 4 Hc cross section using the 

same parameters for the | state at 16.87 MeV with 
E R P = 0.21 MeV as seen in Figure |U However, when 
including modern data as in Descouvemont et al. [52j the 
agreement was found to be poor. Accordingly, we have 
performed a fit of the iS-factor provided by Ref. (HJ that 
gave E* R = 0.35779 MeV, with 7 g c _ 3 = 1.0085 MeV and 
7p = 0.025425 MeV. Hence, our new calculations of the 
3 Hc(d,p) 4 He and 3 H(d,n) 4 He rates coincide with those 
we have been using in previous papers, when the con- 
stants do not vary, showing that the single pole approxi- 
mation is sufficient. 



3 H(d,n) 4 He 



> 




E (MeV) 

FIG. 7: S-factor for the 3 H(d,n) 4 He reaction from Ref. [H 

(dashed) and our calculation using the parameters from 

Barker [k| (solid). See Ref. for the references to ex- 
perimental data. 



D. Reaction rates 

The reaction rates are shown in Figure [9] As expected 
from Eqs. (|4.1j) and (|4.2|) they are only slightly affected 
by variations of Snn (less than 5%). From the sensitiv- 
ity study of Ref. (53, we deduce that 3 H(d,n) 4 He rate 
variations have no effect on BBN while 3 He(d,p) 4 He rate 
variations induce only very small (<4%) changes in the 
7 Li and 3 He abundances. Because the change is so small, 
we can make a linear approximation to the sensitivity as 
shown in Table fl] displaying {8Y/Y)/8mn values for both 
reactions. 

Our results are significantly different from those of 
Berengut et al. [30j : first, we use a more elaborate pa- 
rameterization of the cross-section, second, our calcula- 
tion of the resonance energy shift is less sensitive to the 
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3 He(d,p) 4 He 



.a 




10 ' 3 10 2 10 



E (MeV) 

FIG. 8: As in Fig.[7]but for 3 He(d,p) 4 He and our calculation 
with parameters fit to ref. [52| (dashed). Deviations with 
experimental data at very low energy are due to screening. 



5 NN = -0.03, -0.015, 0., 0.015, 0.030 




10 3 10 2 10 1 1 10 



T(GK) 

FIG. 9: Relative variation of the 3 H(d,n) 4 He (solid line) and 
3 He(d,p) 4 He (dashed line) rates for 5 NN = -0.030, -0.015, 
0.015 and 0.030. 



variation of constants. 

While, for the nominal values of the constants, both 
parameterizations reproduce the experimental data very 
well, their comparison (our Eq. (|4.3[) . and Eq. (9) 
Ref. [30l| ) shows important differences. Only the par- 
tial width energy dependence in the entrance channel 
(Tin = Td) is considered in Ref. [3(J, neglecting the 
outgoing and total energy dependences. To be consis- 
tent, the latter must indeed depend on energy at least 
through T- m . The partial width (T out = T p or T n ) en- 
ergy dependence in the exit channel is small (large out- 
going energy) and can be neglected. But, in addition, 
as shown in Eq. (|4.3p . the total width (r = Td + T p or 
TfZ + r„) energy dependence (essentially from Td) must 



TABLE I: Abundance sensitivity, d logY/<5jvjv, to a variation 
of the N-N interaction at WMAP baryon density. Blank en- 
tries correspond to negligible values. 



Reaction Y p 


D/H 


3 He/H 


7 Li/H 


3 H(d,n) 4 He 






-0.015 


3 He(d,p) 4 He 


-0.027 


-1.14 


-1.10 



be included. We also include in our calculations the en- 
ergy shift AE r (Er). As can be seen in Figure [TUl the 
variation of the S-factor with AEr (now the supplemen- 
tary shift due to variation of constants i.e. ^ AE r (Er)) 
is v ery different from ours when following Berengut et 
al. [30| prescription. (To emphasize the different be- 
haviour, we have explored, in that figure a wider range 
of AE R values ± 135 keV.) 



1 He(d,p) 4 He 




10 2 10 ' l 

E (MeV) 



FIG. 10: S-factor for the 3 He(d,p) 4 He reaction and <5jvat = 
—0.30, —0.15, 0., 0.15 and 0.30 according to our calculation 
(solid black) or following the parameterization of Ref. [3(|] 
with the same energy shifts (dashed pink). To emphasize the 
effect we use larger variations, by a factor of 10, than in Fig. [9] 

These large differences reflect the differences in the pa- 
rameterization of the cross section. This shows that the 
expression of the cross section provided by Cyburt [56| 
and used by Berengut et al. [30l | is excellent for a param- 
eterization of the experimental data, it is by no means 
adapted to a study of the influence of nuclear parameters 
on the cross-section. 

The calculation of AEr as a function of Snn is 
obtained by the difference between the energy of the 
| states and of the two clusters emission thresholds 
in the entrance channel, both depending on the N-N- 
interaction. Berengut et al. [3(| assume that these levels 
follow the dependence of the 5 Li and 5 He ground states 
but these | + state have indeed a 3 Heg)d or tigid structure 
so that our model is more appropriate (see Ref. (57|). 
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Furthermore, we note that significant changes in ele- 
ment abundance's occur when the variation in the quark 
mass, 8m q /m q (their parameter of variation), is of or- 
der 1 %. We recall that Sm q /m q can be related to a 
change in a Yukawa coupling through (1 + S)Sh/h. Us- 
ing Eqs. ([1.5)1 and (|2.5|) we can further relate 5m q /m q 
to Snn (ignoring the contribution from RSa/a) as 

5m q /m q re -OASnn- (4.4) 

Thus for Snn < 0.01, we are essentially restricting 
Sm q /m q < 0.004 and would expect a significantly smaller 
effect on the abundances than seen in Rcf. 1301. 



the Spite plateau is an indication that depletion may not 
have been very effective. 

Astronomical observations of these metal poor halo 
stars [6r| have led to a relative primordial abundance 
of 

Li/H = (1.23±S:?S) x 10~ 10 . (5.4) 

A more recent analysis by Sbordone et al. (67j gives 

Li/H = (1.58 ±0.31) x lO" 10 . (5.5) 

For a recent review of the latest Li observations and their 
different astrophysical aspects, see Ref. (68j . 



V. BBN CONSTRAINTS 
A. Observational constraints 

Deuterium, a very fragile isotope, is systematically de- 
stroyed after BBN. Its most primitive abundance is de- 
termined from the observation of clouds at high redshift, 
along the line of sight of distant quasars. Very few obser- 
vations of these cosmological clouds are available [EH and 
a weighted mean of this data yields a D/H abundance of 

D/H = (3.02 ± 0.23) x 10~ 5 . (5.1) 

In contrast, after BBN, 4 He is produced by stars. Its 
primordial abundance is deduced from observations in 
Hll (ionized hydrogen) regions of low metallicity blue 
compact galaxies. The primordial 4 He abundance Y p 
(mass fraction) is given by the extrapolation to zero 
metallicity but is affected by systematic uncertainties 
such as plasma temperature and underlying stellar ab- 
sorption [5^]. Using the Markov Chain- Monte Carlo 
methods described in Aver et al. |60| and data compiled 
in Izotov et al. [6~D ]. Aver et al. |62l| found 

Y p = 0.2534 ±0.0083. (5.2) 

Given the uncertainty, this value is consistent with the 
BBN prediction. 

3 He on the other hand, is both produced and destroyed 
in stars so that the evolution of its abundance as a func- 
tion is subject to large uncertainties and has only been 
observed in our Galaxy [f53j], 

3 He /H = (1.1 ±0.2) x 10~ 5 . (5.3) 

Consequently, the baryometric status of 3 He is not firmly 
established [641 ] . 

The primordial lithium abundance is deduced from 
observations of low metallicity stars in the halo of our 
Galaxy where the lithium abundance is almost indepen- 
dent of metallicity, displaying a plateau, the so-called 
Spite plateau [65j. This interpretation assumes that 
lithium has not been depleted on the surface of these 
stars, so that the presently observed abundance is pre- 
sumably primordial. The small scatter of values around 



B. Revisiting SBBN with variable couplings 

The results of the former sections can be implemented 
in a BBN code in order to compute the primordial abun- 
dances of the light elements as a function of Snn- We 
can rephrase our analysis of Ref. [22j in terms of 5n n by 
using Eqs. (|1.5[) and (|2.5[) which yields the constraints 

- 0.7% < Snn < +0.5%, (5.6) 

assuming R = 36 and S = 240. To test the importance 
of the variations in the A = 5 rates, we first include only 
those variations in 3 He(d,p) 4 He and 3 H(d,n) 4 He. Fig- 
ure [TT] compares the BBN predictions for different values 
of Snn up to 30%. We emphasize that a 30% variation in 
Snn is unrealistic since it corresponds to a 175% varia- 
tion on Bo- As one can see, the curves for 4 He and D/H 
are nearly horizontal, and the effect on 7 Li is insufficient 
to solve the lithium problem. 

Next, we allow all rates depending on Snn to vary. 
Thus, Figure [12] updates Figure 4 of Ref. [22| assuming 
again S = 240 and R = 36. In this case, we see that 
the D and He abundances can provide constraints on 
the variations of h, compatible with zero, while the 7 Li 
abundance can only be reconciled with observations for 
Sh/h re +3 x 10~ 5 . We now find 

- 2 x 10~ 6 < — < 8 x 10~ 6 . (5.7) 
h 

Recalling that with S = 240 and R = 36, one has 

c 7 

Snn ~ -321 — , (5.8) 
h 

which gives 

- 0.0025 < Snn < 0.0006. (5.9) 

Thus, we see that variations in Snn as large as that 
needed to reconcile 7 Li induce an excess of D/H and a 
deficit of 4 He (even when the large uncertainties in Y p 
are taken into account). 

It is interesting to note that the latter problem can be 
reconciled by adding a relativistic degree of freedom. In- 
deed, it has recently been pointed out that extra relativis- 
tic degrees of freedom may be required from the analysis 
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FIG. 11: Effect of the variation of the N-N interaction induced 
solely by the modification of the nuclear rates of 3 He(d,p) 4 He 
and 3 H(d,n) 4 He on the primordial abundances of the light 
element compared to the constraints obtained in Ref. 



FIG. 12: Update Figure 4 of Ref. [Zj assuming S = 240 and 
R = 36 (solid blue line), using new rates for 3 He(a,7) 7 Li 
[SH and 1 H(n, 7 )D [U and the Q b value from WMAP7 Q. 
The top axis is —Snn from Eq. (|5.8[) (mind the sign) and the 
dashed red line assumes N v = 4. 



of CMB data (zH and BBN data 0. As a consequence, 
we also show in Figure [T2l the resulting abundances when 
number of relativistic species is equivalent to four neu- 
trino families. As one can see, even at large S^n, the 
helium abundance is concordant with observations. In- 
terestingly, although the constraint from 4 He on Sh/h 
for positive values is relaxed, the constraint from D/H is 
tightened and becomes more severe for N v — 4. There- 
fore, to reconcile D and 7 Li abundances we must rely 
on subsequent destruction of D/H to match the quasar 
absorption system data. In fact, high D/H is a general 
consequence of lowering 7 Li from SBBN values [73|]. For 
Sh/h < 0, while the constraint from D/H is relaxed, the 
constraint from 4 He is strengthened and a similar limit 
is obtained. 



VI. DISCUSSION 

This article investigated the influence of the variation 
of the fundamental constants on thepredictions of BBN 
and extended our previous analysis [221 ] . Through our de- 



tailed modeling of the cross-sections we have shown that 
although the variation of the nucleon-nucleon potential 
can greatly affect the triple-a process in stars, its effect 
on BBN and the production of heavier elements such as 
CNO is minimal at best. At the temperatures, densi- 
ties and timescales associated with BBN, the changes in 
the 4 He(aa, 7) 12 C and 8 Be(a,7) 12 C reaction rates are 
not sufficient to produce more than C/H ~ 10~ 21 , and 
is therefore typically 6 orders of magnitude smaller than 
standard model abundances. This conclusion holds even 
when including the possibility that 8 Be can be bound. 
Using the variation of the fundamental constants pro- 
vides a physically motivated and well-defined model to 
allow for stable 8 Be. 

We have also extended previous analysis by including 
effects involving 5 He and 5 Li. This allowed us to revisit 
the constraints obtained in Ref. and in particular to 
show that the effect on the cross-sections remain small 
compared to the induced variation of Bp. This analysis 
demonstrates the robustness of our previous analysis and 
places the understanding of the effect of a variation of 
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fundamental constants on BBN on a safer ground. 

Our analysis can be compared with Ref. [30| who 
reached the conclusion that such variations may increase 
the variation of 7 Li and exacerbate the lithium problem. 
While formally correct, our results show that this re- 
quired large variation of Snn is incom pat ible with the 
BBN constraints. Note also that Ref. 74J assumed an 
independent variation of the energies of the resonances 
while our work considers the variation of the energies of 
these resonance that arises from the same physical origin, 
so that their amplitudes are correlated. 

Finally, we have extended our analysis to include the 
possibility of an extra relativistic degree of freedom. 
Because of the different dependencies, Y(Snn 7 N v ) and 



D/H(8nn, N v ), the limits on 5nn or Sh/h arc not re- 
laxed for 7V„ = 4. The only possibility to reconcile 7 Li 
in this context is a variation of Snn ~ —0.01, along with 
the post BBN destruction of D/H. 
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